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<N ! Abstract 

As a step toward performing a complete coupled-channels analysis of the world data of 
ttN,j*N — > 7rN,rjN,TnrN reactions, the ttN — > 7T7T./V reactions are investigated starting with the 
dynamical coupled-channels model developed in Phys. Rev. C76, 065201 (2007). The channels in- 
cluded are irN, rjN, and nirN which has 7rA, pN, and aN resonant components. The non-resonant 
amplitudes are generated from solving a set of coupled-channels equations with the meson-baryon 
potentials defined by effective Lagrangians. The resonant amplitudes are generated from 16 bare 
excited nucleon (N*) states which are dressed by the non-resonant interactions as constrained by 
the unitarity condition. The data of total cross sections and irN and irir invariant mass distribu- 
tions of Ti + p — ► 7r + 7r + n, 7T + 7r°p and n~p — ► 7r + 7r _ n, 7r _ 7r p, 7r°7T n reactions from threshold to the 
invariant mass W = 2 GeV can be described to a very large extent. We show the importance of 
the coupled-channels effects and the strong interference between the contributions from the ttA, 
aN, and pN channels. The large interference between the resonant and non-resonant amplitudes 
£T^ ' is also demonstrated. Possible future developments are discussed. 
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I. INTRODUCTION 

It has been well recognized |l| that a coupled-channels analysis of the data of meson 
production from irN, jN, and N(e, e') reactions is needed to extract the parameters of 
the excited nucleon (N*) states in the energy region above the A (1232) resonance. This 
has been pursued by using the K-matrix models J2jj [a, M, | a |, the Carne gie- Mellon-Berkeley 



(CMB) model Q and the dynamical models SliMQ [H, M,M, GJ, E]- Since 
two-pion production processes account for about half of the total cross sections of ttN 
and 7./V" reactions, the N* parameters extracted from these analyses are reliable only when 
the employed models are consistent with the two-pion production data. As a step toward 
performing a complete analysis of the world data of ttN, ^*N — > ttN, rjN, ttttN reactions, 
we investigate in this paper the ttN — > irnN reaction starting with the dynamical coupled- 
channels model developed in Ref. [13J (JLMS) from fitting the 11N elastic scattering data. 
A number of theoretical investi gati ons of 7rN — > 7T7tA reactions have been performed using 



(T) tree-diagram calculations [17|, [l8|, [l9|, 120, |2l|, |22j, |23j , (2) chiral perturbation theory [2 



25l . l26l |27| |. an d (3)tree-diagram calculations including coupled- channels effects on nucleon 
To see the scope of our investigation, it is useful to give a brief review of 



resonances 



these previous works in the following subsections. 



A. Tree-diagram calculations 



The tree-diagram calculations are based on phenomenological Lagrangians and the Breit- 
Wigner form parametrization of nucleon resonances. All 7r ± p — > tttcN channels on the 
proton target were investigated in Refs. 2ll . |22J |23j, while only the n~p — > n + n~n channel 
was studied in Refs. 
invariant mass W = 



17, m m 20 



These investigations covered the energy region up to 
1551 MeV (T w = 600 MeV) and only investigated the role of the A(1232) 
and A* (1440) resonances. They could describe reasonably well the experimental data in 
the considered energy region. Attempt was also made to investigate the scalar-isoscalar 
two-pion decay process of A*(1440) — » A(7T7r)g^ ave of the Roper resonance. Although the 
tree-diagram calculation is a convenient tool to catch the qualitative features of the reaction 
processes, it is of course not consistent with the unitarity requirements. Furthermore, such 
an approach starts to break down in the W ^ 1.5 GeV region where more N* states needed 
to be considered. 



B. Chiral perturbation theory 



Chiral perturbation theory calculations of irN — > tcttN reactions have been performed 
up to the order 0(q 3 ). The focus of Refs. [24| was on evaluating the nN — >■ nnN threshold 
amplitudes D\^ and extracting tttt scattering length within the heavy baryon chiral per- 



26| to also compare with the cross 



turbation theory. This work was extended in Refs. [2| 
section data up to W = 1.38 GeV (T n = 400 MeV). They found that the loop corrections 
are small and it is difficult to extract the isoscalar S-wave tttt scattering length because of 
the large uncertainty associated with the iV*(1440) — > A(7T7r)g^ ) ave decay. Similar chiral 
perturbation theory calculation was also performed to compare with more data near the 
threshold by Mobed et al Q. 

Of course, the application of chiral perturbation theory to investigate the ttN —>■ inrN 



reactions near threshold is an important advance. But it is not clear how it can be used to 
investigate the nucleon resonances up to W = 2 GeV. 



Tree-diagram calculation including coupled-channels effects on N* 



The tiN — > irirN study performed by the Julich group [28| also used the tree-diagrams 
generated from phenomenological Lagrangians, but including the P33(1232), Pn(1440), 
.Di3(1520), Sii(1535), and S f n(1650) resonance states. They focused on the low W < 1.38 
GeV (T n = 400 MeV) region. The coupled-channels effects are included by using the non- 
resonant amplitudes generated from the ttN coupled-channels model of Ref. [8J to evaluate 
the self-energy and the decay functions of the considered resonances. They could describe 
to a large extent most of the available data in the W < 1.38 GeV region. They found that 
their calculations started to break down in the higher W region. 



In this work, we depart from all of the earlier works described above by considering the 
whole energy region from threshold up to W = 2 GeV and all nucleon resonances listed by 
Particle Data Group [34j (PDG). The calculations are performed by using the JLMS model. 
Schematically, the following coupled integral equations in each partial wave are solved within 
the JLMS model 

T a ,p(p a ,Pfi;E) = VafiipcPp) + J2 P 2d P V a,s(p a ,P)Gs(p,E)T S)f3 (p,p l3 ,E) , (1) 

s J 

with 

Va,0{p a ,P0) = VajipnPp) + ^ ~^ TJq , (2) 

N* ^ ~ lV1 N* 

where a, j3,5 = jN, tcN, rjN, and titcN which has 7rA, pN, aN resonant components, Gs(p, E) 
is the propagator of channel S, M^* is the mass of a bare excited nucleon state iV*, v a $ is 
defined by meson-exchange mechanisms, and the vertex interaction Tn*,p defines the N* -^ f3 
decay. The fits to irN elastic scattering data were achieved by including one or two bare iV* 
states in all S, P, D and F partial waves. The details can be found in Ref. [l3j. 

While the non-resonant interactions v at p are deduced from phenomenological Lagrangians, 
the model contains many parameters mainly due to the lack of sound theoretical guidance 
in parametrizing the bare iV* — > ttN, rjN, 7r A, pN, aN form factors. Although fitting the tcN 



elastic scattering data is rather complex and time consuming, as reported in Ref. [131 ] . it is 
clearly not sufficient to determine these parameters; in particular the parameters associated 
with the unstable 7rA, pN, and aN channels. Thus it is important to test the JLMS model 
in the study of 7riV — > inrN reactions which are known to be dominated by these channels. 
This is the main purpose of this work. The results presented in this paper are therefore 
obtained using the parameters taken from the JLMS model. 

It is tempting to try to use the available 7riV — > irnN data to improve the JLMS model. 
However this is not an easy task at the present time mainly because of the lack of sufficient 
experimental data in the considered energy region up to W = 2 GeV (T n ~ 1.6) GeV. 
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FIG. 1: Graphical representations of T^n^-kN of Eqs.([3|)-([7|). 

Furthermore most of the data at high W, obtained before 1970's, are not accessible [35(. In 
an effort by R. Arndt [35], some of the final irn and irN invariant mass distribution data 
have been recovered and considered in this paper along with the most studied total cross 
section data [33j. This allows us to investigate these two observables for all possible final 
7nrN states of n p reactions and in the entire energy region from the threshold to W = 2 
GeV. Very limited data on the angular distributions can be found in the literature, only few 
in the low W region and practically nothing in the higher W region where we need pin down 
the parameters associated with many iV* states. We thus will not consider these observables. 
For the same reason it is difficult to use the available irN — > ttttN data at present time. 
This is mainly due to the problem that the minimization in determining a large amount 
of parameters of our model is heavily weighted by the very extensive and far more precise 
data of iiN elastic scattering. We thus focus in this paper on investigating the dynamical 
content of the JLMS model, in particular on several aspects which were not addressed before 
such as the effects of coupled-channels and the role played by the interference between the 
different channels. Even with this limitation our investigation is clearly more extensive than 
all previous works reviewed above. 

In section II, we briefly recall the formulas of Ref. 12) for calculating the ttN — > mrN 



amplitudes within the JLMS model. To give more details about our calculations, explicit 
expressions for calculating the total cross sections and tcN and tttt invariant mass distribu- 
tions are given in section III. The results are presented in section IV. In section V, we give 
a summary and discuss possible future developments. 



II. FORMULATION 



Within a Hamiltonian formulation [12J within which the JLMS model was developed, the 
ttN — > ttttN amplitude is illustrated in Figfj] and can be written as 



with 



T- kN ,-kttn{.E) - T^^ N {E) + T* N ^ N {E) + T£ Nj7T7tN (E) + T% N ^ wN (E), (3) 

T n N, n7T N(E) = v ttN,ttitN + 2 , Tt T n,Mb{E)Gmb{E)vmB,ttttN , (4) 



MB 



TkN,-k-kn{E) = T n N tn A(E)G 7r ^.(E)T A-*nN , (5) 

TwN,ttttn( e ) — T 7rNtpN (E)G pN (E)h p ^ 7T7r , (6) 

TnN,imN\E) = T nNt(rN (E)G (T N(E)h (T _ tn7T , (7) 

where Ta^ttAt, h p ^ n7r , and h^^ describe the A(1232) — > ttN, p — > 7T7t, and a — »■ 7r7r decays, 
respectively, Gmb{E) for MS = 7iN,r]N,TcA, pN,aN are the meson-baryon propagators. 
The 7riV — > Mi? transition amplitudes are 

T-kn,mb{E) = t nN ^ MB (E) +t 7TNMB (E) . (8) 

The second term in the right-hand-side of Eq.pl) is the resonant term defined by 

t MB,M'B l( y E ) = Yl ^MB-^N*{E)[D(E)] i} jT N *^M'B'{E), (9) 

JV* jy* 

with 

\D{E)-\ 3 {E) = (E- M° N: )5 hl - ttj(B) , (10) 

where Mjy* is the bare mass of the excited nucleon state N*, and the self-energies are 

Sjj(-E) = Y Fn*^mbGmb{E)T mb ^ n *(E) . (11) 

MB 

The dressed vertex interactions in Eq.flU]) and Eq.flTTJ are (defining Tmb-*n* — ^n*-*mb) 

FmB->N*(E) = T MB ^ N * + 2^ tMB,M'B'{E)GM'B'{E)T M i B i^ N * , (12) 

M'B 1 

Fn*-^Mb{E) = Tn*^A1B + 2^ T N*^M> B'G m> B'{E)tM> B> ,Mb{E) . (13) 

Af'B' 

The non-resonant amplitudes tMB,M'B' in Eq.® and Eqs.(IT2"l)-(IT3"l) are defined by the fol- 
lowing coupled-channels equations 

tMB,M',B'(E) = VMB,M'B'{E) + 2_, v MB,M"B"(E)G M "B"(E)tM"B",M'B'(E). (14) 

Af"S" 

The channels included are MS = 7riV, 77JV, 7rA, pN, aN. All quantities defined above are 
described in detail in Refs. [12J, ll3j and can be calculated within the JLMS model. The 



only exception is the direct production term v^n^n in Eq.flU) which is not included in the 
JLMS model. The procedure for deriving ty y^ /y from Lagrangians by using the method of 
unitary transformation is explained in Ref. 12J. In this work we consider a v^n^-kN model 



which involves only the N and A(1232) intermediate baryon states such that its parameters 
can also be consistently taken from the JLMS model. Our calculations thus do not have any 
adjustable parameters. The mechanisms of the considered v n N,mrN are illustrated in FigfJl 
The calculations of the terms T^ N B ^ N with MB = nA,pN,aN, defined by Eqs.©-©, 
are straightforward. On the other hand, the calculation of the second term of T^§ nnN , 
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FIG. 2: The considered v w 



N,ttttN- 



defined by Eg. (HI), is much more complex. To simplify the calculation, we first note that the 
mechanisms (a)-(e) in the upper row of Fig. [2] can be written as 

V ttN,-kttN ~ V 7T N,nNG n N(E)hN^TTN, (15) 

where v^ N ^ N is the sum of all mechanisms (a) - (e), v^n^n is part of the driving term in 
calculating the amplitude TnN,irN{E) from JLMS model, and tiN^nN is the N — > ttN vertex 
function. If we neglect the final state interactions on the mechanisms (f)-(k) of FigfSJ we 
can write 



^Jv.TTTrTV ~ V ttN.tt7tN + \ v nN.nN + 2 .T 7T N,MBGMB(E)vMB,-irN{E)]G 7T N(E)h N ^ 7T 



N 



MB 



,(/-*) 



~ v { J N ^' nN + T nN ^ N (E)G nN (E)h N ^ nN . (16) 

In deriving the above equation, we have used the full scattering equation T = V + TGV 
of the JLMS model with the approximation that only the A (1232) bare state and the ttN 
state are kept in summing the intermediate states. We use Eq.( TT6l) in this work. 



III. CROSS SECTION FORMULA 

For presenting results, we here give explicit formula for calculating the cross sections of 



ttN — > ttttN reactions. Within the formulation of Refs. jl2|, [l_3[ , the S-matrix is defined by 

S fl = 6 fi - 2ixi5\P f - Pi)(f\T\i) , (17) 

and the plane-wave state is normalized as 

(q\p) = 5(q-p). (18) 

In the center of mass frame, the momentum variables of the 7iN — » ttttN reaction with 
invariant mass W can be specified as 

a(Pa) + KPb) -> c{Pc) + d(p d ) + e(p e ) (19) 

where p a = —pb = k with k defined by W = E a {k) + E b (k) , p c +Pd = —p e = k', and (c+d + e) 
can be any possible charged states formed from two pions and one nucleon. The total cross 



section of the process Eq. (TT9|) can then be written as 

Oab^cde = -(2vr) 4 / dp c dp d dp e 5 4 (p a + Pb ~ Pc ~ Pd ~ Pe) 

V J 

J2\(PcPdPeJ\ T \PaPb,i}\' 



(2s a + l)(2s b + l)f- F 
E a (k)E h (k) 4 r dp c dp d dp e 4 

wk (27r) Jx^^ 6{Pa + Pb ~ Pc ~ Pd ~ Pe 

1 



x (2ga + 1)(2gfe + 1) E I V 'E c E d E e (p c PdPe, f\T\p a pb, i)\ 2 (20) 

where i, f denote all spin (s a ,s az ) and isospin (t a ,t az ) quantum numbers, and Ykj means 
summing over only spin quantum numbers. The above equation can be written as 

rW—rrie d(J 
Oab^cde = / Ul~ dMcd ( 21 ) 



with 



dM, 



C(l 



(22) 



Here we have defined 



kE a {k)E b {k) 

Pi — 7T , 

W = E e {k') + E cd (k') 



E cd (k') = ^Ml d + ik'f , (23) 

M cd = E c (k cd ) + E d (k cd ) , 

and k cd is the relative momentum between c and d in the e rest frame. 

The T-matrix elements in the above equation are calculated from Eqs. (j4"]) -f7j). For the 
quasi two-body processes Eqs.(E])-(ED w hh a resonant unstable particle R which decays into 
c + d state, the T-matrix elements take the following form 

, ,, ,, f\T\h '\ V^ (Pc s cz , t cz ; p d , Sdz, t dz \Hj\k', s Rz , t Rz ) 

[p c p dPe ,J\i\K,t) - ^ E-E e (k')-E R (k')-E eR (k',W) 

x \k , s Rz , t Rz ] k , s ez , t ez \l I /c, s az , t az ] k, sj) Z , ifoz/- (24 J 

For any spins and isospins and cm. momenta p and p', the MB — > M'B' T-matrix elements 
are in general defined by 

(P ,SM'z,tM'z] ~P > S B > z ,t B 'z\T\p, SMzi^Mz] ~P, S Bz ,tB z ) 

= E Yl>,m> l (P')YIm l (p) 

S,S',L,L',J,M,I 

X(sjVf') %'z, Sfi/, SB'^|5' , S Z )(L , M L , S 1 , SJ J, M)(tM',tM' Z ,tB',tB> Z \I, h) 
x (sm, %2, Sfi, Ss^S", S Z )(L, Ml, S, S z \J, M)(tM, tjtfz, ^B, tBz\I, I z ) 

xlv',S'L'J\T\W)\p,SLJ), (25) 



where (ji,ji,mi,m2\jm) is the Clebsch-Gordan coefficient of the j\ + ji — J coupling, 
(k f , S'L , J\T I (W)\k, SLJ) is the partial wave amplitude defined by the total angular mo- 
mentum J, orbital angular momentum L, total spin S, and total isospin /. We use the 
JLMS model to generate (k', S"L'J|r J (W)|A;, SLJ). 
As for the R —>■ cd vertex function, the combination 



E c (p c )E d (p d )E R (k')(p c , s cz , t cz ; p d , s dz , t dz \Hj\k', s Rz , t Rz ) 



(26) 



is Lorentz invariant and therefore can be written in terms of its matrix element in the rest 
frame of R 



(p c , s cz ,t cz ;p d , s dz , t dz \Hj\k , s Rz ,t Rz ) 
S(p c + Pd-k' 



E c (k cd )E d (k cd )M R 



\ E c (p c )E d (p d )E R (k 



y:{k c d, Scz, t cz ] —k c d, s dz , t dz \Hi\0, s Rz ,t Rz ), (27) 



with 



ci s czi s di s dz\S c di S cdz ) [L cd , 7Tl cd , S cd , S cdz \S R , S Rz ) 



\kcdi s czi tcz] &cdi s dzj t dz \lij\0, S Rz: t Rz ) 

= E K* 

L c d, m cd,Scd,S c dz 

x (t c , tcz, td, tdz\tR, ^fi2)5 / L cd ,m cd (^cd)-^I,« ,S« {k c d)}- 

The vertex functions are 

TT /V I TT /V 






where L$ N 



1 9 A 



3/2, L' m 



0,^ 



0, LL 



1 S p 

^1 U TTTT 



(28) 

(29) 

(30) 

(31) 

1. Here it is noted that the 



factor v2 in Eqs. (l3"0]) - (f3~Tj) comes from the Bose symmetry of pions, and the phase factor 
% and (—1) are chosen to be consistent with the non-resonant interactions involving tcNA, 
aim and pmr vertex interactions defined in Ref. 12J. The form factors in Eqs. (l30l) - (I31I) 
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m* (2vr)3/2V 3^^)^ 2E N (q) 
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fjo 
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with f nNA = 2.049, g 



0.7750, g^ = 0.6684, A wNA = 649 MeV, A 



, (32) 
(33) 

(34) 
378 MeV, 



and A 



pmr 



461 MeV. The above vertex functions are determined from fitting the nN phase 



shifts in P33 29J] and mr phase shifts [30j. 

With the vertex function f R ^ cd (q) given above, the self-energy appearing in eR Green 
function S e ^ (eR = irA,Np,Ncr), are calculated from (see Appendix IA1 

M eR (q) \fR->cd(q)\ 2 



Z eR (k,E) 



m R 
E R (k) 



q dq 



[M* R (q) + F]V2 E _ Ee (k) - {[E c (q) + E d (q)f + fc 2 }V2 + lE 

(35) 
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FIG. 3: The total cross sections predicted (solid curves) from the JLMS model are compared with 

7rN,7T7vN' 



the data [33]. The dotted curves are from turning off the amplitude T 



To derive the above equation, we have used the Lorentz transformation to calculate the 
self-energy in arbitrary frame from the vertex function defined in the rest frame of R. 

IV. RESULTS 

With the formula presented in sections II and III, we now present results for the total 
cross sections <7 nN __> wnN and the invariant mass distributions da/dM nN or da/dM n7T for all 
possible final ttttN states of ii ± p reactions in the energy region from threshold to invariant 
mass W = 2 GeV. As mentioned in section I, our investigation thus is much more extensive 
than all of the previous dynamical calculations of ttN — » ttttN reactions in both the energy 
range covered and the N* states considered. 

We can predict the ttN — > ttttN cross sections using the information generated from 
the JLMS model except that due to the direct production interaction v n N t7T7T N in Eq.(j3]). 
This term, which can induce the ttN — > rnrN — > ttN mechanism to influence the ttN 
elastic scattering, was not included in the development of JLMS model. We thus first 
examine the importance of this term. As discussed in section II, the contributions from 
v^n^ttNi calculated by using T^ n7rN defined in Eq.f TPoj) . is completely fixed by using the 
same parameters from JLMS model. Thus no additional parameters are introduced in our 
calculations. 

Our results for the total cross sections are shown in Figj3j The solid curves are from 
our full calculations and the dashed curves are from turning off the term T^ nnN . We 
see that both the magnitudes and the energy-dependence of the data for all five two-pion 
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FIG. 4: The total cross sections predicted (solid curves) from the JLMS model are compared with 
the data [33]. The dotted curves are from turning off the amplitude T^ n7rN . 

production processes can be reproduced to a very large extent by our full calculations (solid 
curves). Clearly, the direct ttN — > tctcN mechanisms play a significant role in determining 
the predicted cross sections. In particular, it is instrumental in obtaining the agreement 
with the n + p — > 7i + 7T + n data. Its effects at low W can be more clearly seen in Fig. HI Here 
we also see that the agreement with the data of n^p — > 7r°7r~p is mainly due to the effects 
of v n N,irirN- We note here that our full calculations (solid curves) in FigH] are comparable 



to those of the chiral perturbation theory calculation of Refs. [23, |26|, |27|]. This suggests 
that the model v n N,mrN considered here is fairly reasonable and the discrepancies with the 
data in the higher W region, seen in FigJ21 are more likely from the uncertainties in the 
contributions from irA, aN, pN transitions. 

The main feature of this investigation is a dynamical coupled-channels treatment of ttttN 
channel which has the irA, pN, aN resonant channels. In our calculation, this effect can be 
explicitly seen by writing the coupled-channels equations, Eq. (TT4T) . as 



t 



nN,MB 



(E) 



22 l l -v^knm'b'^m'b^mb 
m'b' 



v«« 



(36) 



where MB = nA, pN,aN, and the intermediate meson-baryon states can be M B = 
ttN, r]N, ttA, aN, pN. When only the term with M'B' = MB in the Eq.([3S]) and in Y N *^ MB 
of Eqs. (TT2"|) - (TT3"j) is kept, the calculated total cross sections (full curves) are changed to 
the dotted curves in Figf5J If we further neglect the coupled-channels effects by setting 
t?rN,MB = v nN,MB, we then get the dashed curves which are very different from the full 
calculations (solid curves), in particular in the high W region. 

To see the coupled-channels effects more clearly, we show the corresponding results for the 
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States 


W (MeV) 


"total V ' 






ttN 


rjN 


7rA 


oN 


P N 


Sn 


1535 


50.7 


34.9 


8.0 


6.1 


0.3 




1650 


72.3 


4.0 


1.1 


8.2 


14.4 


#31 


1620 


26.6 





67.3 





3.1 


Pll 


1440 


62.5 





7.7 


29.5 


0.3 




1710 


51.2 


2.0 


8.23 


24.6 


14.0 


Pl3 


1720 


27.9 


0.2 


70.2 


0.9 


0.8 


PSI 


1910 


79.5 





6.0 





14.4 


P33 


1232 


100 
















1600 


76.4 





13.6 





10.0 


D 13 


1520 


58.0 


0.0 


36.8 


4.5 


0.8 


D 15 


1675 


38.6 


0.7 


56.5 


0.7 


3.5 


D33 


1700 


12.4 





85.2 





2.4 


*i 5 


1680 


70.0 


0.0 


4.8 


19.3 


5.8 


-F35 


1905 


10.9 





51.1 





38.0 


^37 


1950 


39.8 





59.7 





0.5 



TABLE I: Branching ratios of the 7riV — ► 7riV, rjN, ttA, aN, pN partial wave cross sections calculated 
from the resonant amplitude Eq.©. W is the total energy. 

it N and inr invariant mass distributions at W = 1.79 GeV in FigsJEE Our full calculations 
(solid curves) are able to reproduce the main features of the data. Comparing them with the 
dotted and dashed curves, it is clear that the coupled-channels effects can change strongly 
both the magnitudes and energy- dependence of the nN -^ 7r7riV cross sections. 

To further see the dynamical content of our model, we show in FigsJHJl9]the contributions 
to the invariant mass distributions at W = 1.79 GeV from each of the processes via the final 
MB = 7rA, pN, aN defined by Eqs.Q-©. The results shown in FigsJHHH] indicate that the 
full coupled-channels calculations (solid curves) involve rather complex interference effects 
between these three unstable particle channels. To improve the model, we need to tune their 
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FIG. 5: The coupled-channels effects on ttN — > 7T7T./V reactions. The solid curves are from full 
calculations, the dotted curves are from keeping only M'B 1 = MB in the Eq. ([36|) and in Tm* ^-m p 
of Eqs. (fT2]) - (fT3]) . the dashed curves are from setting tMB,M'B' = vmb,M'B'- The data are from [35]. 



relative importance. 

Compared with all previous investigations, another feature of this investigation is our 
treatment of the N* resonance amplitudes. These amplitudes are generated from 16 bare 
states, as given in Ref. [131], which are dressed by the non-resonant interactions, as required 
by the unitarity condition and denned by Eq.(j9"|)- ([T51) . In Fig. [TO], we compare the full 
results (solid curves) and that calculated from keeping only the non-resonant amplitudes 
(dashed curves) for the invariant mass distributions of ir~p — > TT + 7r~n at W = 1.44, 1.60, 
1.79 GeV. Here we note that the peaks of dashed curves in M(ir + n) and M(ir~n) distribu- 
tions around 1.2 GeV in Fig. [TO] are due to the decay of A in the intermediate 7rA state of 
the "non-resonant processes" of ttN — > ttA — ► ttttN, whose amplitude is defined by Eq. ([5]) 
but replacing T^n-A with its non- resonant amplitude t^^^A generated from Eq. (fT4"l) . Sim- 
ilarly, the decay of p (a) in the intermediate pN (crN) of the "non-resonant processes" of 
ttN — > pN(aN) — > ttttN can be responsible for the peaks of dashed curves in the M(7r + 7r~) 
distributions (lower row of Fig. [TO]) . 

By comparing the solid and dashed curves in Fig. [TO] it is clear that the full calculations 
involve comparable contributions from resonant and non-resonant amplitudes. In the same 
figure, we also show ttttN phase-space distributions (dotted curves) normalized to data. The 
shapes of both theoretical results deviate significantly from the phase-space. 

The results shown in Fig. [TO] indicate that the fits to the data depend strongly on the 
parameters associated with the N* — > 7rA, pN, aN vertex functions which are treated 
purely phenomenologically within JLMS model. One possible improvement of the model 
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FIG. 6: Coupled-channels effects on the invariant mass distributions of n + p — ► 7T7tA at W = 
1.79 GeV. The solid curves are from full calculations, the dotted curves are from keeping only 
M'B' = MB in the Eg. ([36]) and in Tn*^mb of Eqs. ([T2|) - (fT3|) . the dashed curves are from setting 
tMB,M'B' = vmb,M'B'- The data are from R. Arndt [351 ] . 



is to explore how these vertex functions can be calculated from sound hadron structure 
calculations. An attempt along this line for a two-channel ttN scattering in Sn state was 
pursued in Ref. [311] using the constituent quark model, but was not successful. 

The complexity of the calculated resonant amplitudes can be further seen in Table I 
where we show the calculated branching ratios of the contributions from each channel to the 
partial- wave cross sections calculated from the resonant amplitude t^ N MB (Eq.fjHJ)) at the 
resonant energies listed by Particle Data Group. Clearly, the resonant amplitudes involve 
strong interference between the irN — > 7r A, pN, aN — > irnN amplitudes. To improve the 
model, we need to tune their relative importance. Clearly more detailed data, such as 
the single or double angular distributions and polarizations, are needed to make significant 
progress. We emphasize here that the results listed in Table I are not the branching ratios 
of A^* — > MB decay widths at the resonance poles which will be extracted from using the 
analytic continuation methods developed in Ref. [32f] . These results just give some ideas 
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FIG. 7: Coupled-channels effects on the invariant mass distributions of n~p — ► 7T7T./V at TV = 
1.79 GeV. The solid curves are from full calculations, the dotted curves are from keeping only 
M'B' = MB in the Eg. ([36]) and in Tn*^mb of Eqs. (fT2]) - ([T3"]) , the dashed curves are from setting 
tMB,M'B' = vmb,M'B'- The data are from R. Arndt [351 ] . 



about the relative importance between different channels at some energies listed by PDG. 

In Fig. [10], we also observe that our predictions do poorly in describing the tc + -k~ distri- 
bution at low W = 1.44 GeV. We have found that this is the case for all two-pion invariant 
mass distributions of 7i~p — * TT + 7i~n and 7r~p — > Tr°ii n reactions at low W < 1.5 GeV. This 
is given in more detail in Fig. [TTJfor n~p — ► TT°7i°n reaction. We see that our prediction (solid 
curve) does not reproduce the data from the Crystal Ball collaboration [36(. We also show 
the results from keeping only the non-resonance amplitudes (dashed curve) and only the 
resonant amplitude (dot-dashed curve). The shapes of all theoretical curves are similar to 
phase-space (dotted curve) and are far from the data. We have found that the problem can 
not be easily resolved by simply adjusting iV* parameters, in particular those in the most 
controversial P\\ partial waves. It requires detailed analysis and more extensive ttN — > mrN 
data to resolve the problem. 
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FIG. 8: Contributions from ttA (dashed) and pN (dot-dashed) channels to the invariant mass 
distributions of 7r + p — ► nirN at W = 1.79 GeV. The data are from R. Arndt \3l 



V. SUMMARY 



In this paper we have shown that the predictions from the JLMS model can describe 
to a rather large extent the data of total cross sections and nN and tttc invariant mass 
distributions of tcN — > ttttN reactions in the energy region from threshold to W = 2 GeV. 
Our investigation is thus more extensive than all previous dynamical calculations of this 
reaction in both the energy range covered and the N* states considered. 

We have demonstrated the importance of the coupled- channels effects and strong inter- 
ference between the 7rA, pN, and aN. The problem in identifying the mechanisms for 
improving the considered JLMS model is further complicated by the finding that the con- 
tributions from resonant and non-resonant amplitudes are comparable. 

An important finding in this work is that the direct v^n^n mechanisms, illustrated in 
Fig. [21 play a significant role in obtaining good agreement with the data, especially in the 
W < 1.4 GeV region where our results are comparable to those from the chiral perturbation 



theory calculations [25|, |26|, |27J]. This raises the question on the extent to which our elastic 
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FIG. 9: Contributions from ttA (dashed), aN (dotted) and pN (dot-dashed) channels to the 



invariant mass distributions of tv p — ► ttttN at W = 1.79 GeV. The data are from R. Arndt 35]. 



scattering results will be changed by the effect due to v^^irirN an d how the unitarity condition 
is violated. For the former one, it can not be answered easily since it will involve solving 
three-body tcttN — ► tcttN scattering equations, as discussed in Ref. 12|. We however can 
examine the unitarity condition by comparing the total cross sections calculated from (a) 
using the optical theorem Copt = (27r /k)F wN ^ N (6 = 0), (b) er (tot) = a^N^N+a^N^+a^N^N 
(no T^ j7r7r7V ), (c) a (tot) = <T 1 rN,i r N+<r*N,irri+<7irN,m(N- These are shown in FigfTJ Here, a^ N ^ N 
and (J n N,nr] are directly calculated from the employed amplitudes generated from JLMS model 
within which the effects of v n N,mrN are not included, a n N,mrN (no T^§ nvN ) and cr n N,mrN are 
calculated from using the formula given in section II and III, with T^ n7rN defined by v^n^n 
through Eq.(16). We see that (a) (solid curves) and (b) (dotted curves) agree completely 
as required by the unitarity condition within the JLMS model. Their differences with (c) 
(dashed curves) measure the violation of the unitarity condition when the effects due to 
Vkn^kN are n °t consistently included in solving the coupled-channels scattering equations. 
Clearly, the unitarity condition is violated significantly mainly in the high W region. For 
example, the results for n~p total cross sections at W = 1.8 GeV are: (a)=32.94 mb, 
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FIG. 10: The invariant mass distributions of n~p — > 7r + 7r~n at W = 1.44, 1.60, 1.79 GeV. The 
solid curves are the full calculations, the dashed curves are from non-resonant amplitudes, and the 
dotted curves are the phase-space normalized to the data. The data are from Ref. [351 ]. 



(b)=32.48mb, and (c)=35.57 mb. However, the results shown in Fig{12]just mean that 
the effect of v^n^n wm n °t change significantly the elastic differential cross sections at 
forward angles. For a complete unitary calculations for all ttN reaction observables, we 
need to include v n N,mrN effects in solving the coupled-channels equations, as detailed in 
Ref. 12J. This is being pursued along with our effort in developing a combined fit to the 
world data of nN, 'yN — > nN, rjN, ttttN. Our progress in this direction will be reported 
elsewhere. 

Our analysis presented in Figs. loTfTTl indicates the complication of the ttN — > ttttN prob- 
lem. To improve our model, more experimental data, such as the single or double angular 
distributions and polarization observables, are needed to pin down the parameters of the 
model. With the recent effort [35], progress in this direction could be realized in the near 
future. Of course, experimental efforts at the new hadron facilities such as JPARC are highly 
desirable. 
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FIG. 11: The invariant mass distributions of ir~p — ► 7r°7r°n at W = 1.40 GeV. The solid curve is the 
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are from Ref. [36|. (The data is transformed to da/dM using the relation da/dM = 2Mda/dM 2 .) 
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APPENDIX A: SELF ENERGY IN UNSTABLE MB PROPAGATORS 



In this appendix we give a derivation of Eq.(35). To be more explicit, let us consider 
eR = 7rA for Eq.(35) and suppress spin-isospin indices. The starting point is Eq.(21) of 
Ref. 12( which defines the formulation used in JLMS model and this work. 

Since the vertex interaction Hj = T^at conserves the total three momentum of the 
system, we have 



$aP'JHj 



P-k-kN 



-HiIpaPit) 



un-connected 



E - K„ - K v - K N + ie 

= 5(p' A - PA)S(p' w - &)E,rA(Pir, E). (Al) 

The kinematics for evaluating Eq.(Al) is illustrated in FigJTB"! To proceed further, we then 
use the following well known relativistic kinematic relations (for example, see Ref. 38] and 
section 2.3 of Ref. 0) 



P 

Q 



rCjr ~\~ tC 



A'- 



(A2) 



2M nN (E 7rN + M„ 



N, 



x \(Ml N + 2E N (k N )M nN + m 2 N - m 2 n )k n - {M 2 N + 2E n (K)M vN + ml - m 2 N )k 



n 



with 



We then have 



with 



EnN = E 7T (k 7r ) + E N (k N ), 


M nN = E n (q)+E N (q) 


= [El N -P 2 ] l/2 - 


— * — * 

dk n dkN = 


d(k„, k N ) 
d(P, q) 


dPdq, 


d(k v ,k N ) En{K)E N (k N ) 


EM + E N {q) 
) E w (q)E N (q) 


d{P,q) E n (k w )+E N (k N 


E n {k 7T/ 


E N (k N ) 


M wN (q) 



E^qjE^q) [M 2 N (q)+P^ 



(A3) 



(A4) 
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(A6) 



"V 




FIG. 13: Graphical illustration of Eq.(Al) for calculating the A self energy in ttA propagator. 
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With the Lorentz invariance property Eq.(27), we can calculate the matrix element of Hi in 
terms of that in the rest frame of A 



(Mwv|#j|pa) = $(pa -K-k 



N, 



\ E 7r (k n )E N (k N )E A (p A ) 



By using the above relations, we then have 

P 

{p'aP'JHi- ™ N - — —H^p^) 

E — K n — K n — K N + xe 

= 6(p w -p„) (p'^H^kN)- — — n " — r r&*k N \Hi\p A ). 

J E - E n (p n ) - E n (k n ) - E N (k N ) + it 

(A* 

By using Eqs.(A5)-(A7), we then obtain 

P 

(p' A P\\Hi- ™ N — —H^p^) 

E — K n — K n — K N + xe 

- '<*. -a) I} tS&^ w -«i*ib> W. - m>. - P) 

{ E n {K)E N {k N ) M nN (q) 1 dqdP 

X 



E,(f)E K ($) [Jl^ w («) + P 2 ] 1/2 £ - £„(A) - {[£„® + -Bn(?)] 2 + P 2 } 1/2 + it 

, y M„„(g) dm-<J\Hi\0)\ 2 (Aq) 

J [Ml N {q)+P*y/*E-E^)-{[EM+E N m 2 +piy/ 2 + ie- l ] 

Comparing Eqs.(A9) and (Al) and using the partial wave expansion Eq.(28), we then obtain 
in the center of mass frame (p n = —p A = k) 

y (h F) = mA f 2 d M ^(q) I/a-^vI 2 

nA{ ' ' E A (k) J q q [M 2 N (q) + k 2 }V 2 E - E^k) - {[EM + E N {q)Y + k 2 Y' 2 + ie 

(A10) 

Eq.(AlO) is Eq.(35) for eR = ttA. 
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